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Abstract. Scintillation observations of the nearby single millisecond pulsar (MSP) PSR J0030+0451 were carried 
out with the Parkes 64m radiotelescope at three different epochs in 1999. From analysis of the dynamic spectrum 
we obtained the amplitude of the electron density power spectrum logC 2 ^ ~ —3.33 and a scintillation velocity 
Viss Ss 15 km s _1 . This result shows that the Shklovskii effect on the spin-down rate P is negligible. We also 
performed a correlation analysis between pulsar proper motions (Vpm) and scintillation velocities (Viss) using 
updated measurements for a sample of 77 objects, 17 of which are MSPs. The full sample shows a correlation 
coefficient r 3 ~ 80% at an extremely high significance level, while for the MSP sub-sample (excluding 2 outliers) 
we obtain r s ~ 90%. 
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1. Introduction 

PSR J0030+0451 was ind ependently d iscovered by the 
\ Arecibo Drift Scan Search ( Bomcr 2000|) and the Bologna 
sub- millisecond pulsar survey ( D'Amico 2000 ). It was 
' the first millisecond pulsar (MSP) discovered with the 
Northern Cross radiotelescope and it appeared to be 
strongly scintillating at A = 70 cm. Follow-up observations 
performed at Parkes with the 21 cm multibeam receiver 
confirmed the pulse period P ~ 4.865 ms, the dispersion 
measure (DM) ~ 4.33 pc cm~ 3 and showed the pulsar is 
not in a binary system. Updated ephemeris from the on- 
going monitoring were used for the data analysis of our 
scintillation observations. 

Given its galactic coordinates I ~ 113°, b ~ —57° 
and DM, the Taylor & Cordes (1993) electron distribution 
model gives a distance D = 0.230 kpc. Timing observa- 
tions of PSR J0030+0451were carried out by Lommen et 
al. (2000) using the Arecibo radiotelescope. Among the 
derived parameters, they report P ~ (1.0 ± 0.2) x 10 -20 
s s~ x , and an upper limit on the pulsar proper motion of 
H ^ (60 4- 70) mas yr^ 1 . The so called "Shklovskii effect" 
could affect the intrinsic P by 



where V± is the pulsar transverse velocity in km 
s _1 and Z?kpc its distance in kpc ( |5hklovskii 1970 ; 
pamilo et al. 199*4 ). This puts an upper limit to the 
proper motion of /i 60 mas yr -1 . 

MSPs have an observed mean transverse veloc- 
ity of 85 ± 13 km s _1 , a factor ~ 4 lower than 
that of ordinary pulsars, with those in binary sys- 
tems having about twice the velocity of the iso- 



latcd ones (Toscano et al. 1999b Johnston et al. 1998 



Nicastro fc Johnston 1995| b Improving this statistical 
sample is particularly important in the light of neutron 
star birth scenarios. PSR J0030+0451's distance, its pe- 
culia r position in the Loc al Interstellar Medium (LISM) 
(see Toscano et al. 1999a ) and its strong flux variations 
make this MSP a good target for scintillation studies. 

The correlation between pulsar proper motions (V pm ) 
and scintillation speeds (Vj ss ) was investigated in the 
past by several authors (e.g. Harrison & Lyne 1993; 
pupta 1995 ) showing the latter is a good estimate of the 



Pa/P 



1 „-l 



1.1 x 10- 22 VlD kpc s 



(1) 



pulsar transverse velocity. Here we revisit the analysis on 
a sample of 77 objects, including PSR J0030+0451 and 16 
more MSPs, for which proper motion (13 are high confi- 
dence upper limits) and scintillation measuremnts exist. 
The analysis was performed both using the whole sample 
and the MSPs subsample only. 
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Fig. 1. Dynamic spectrum of PSR J0030+0451 with a frequency resolution of 0.5 MHz and a subintegration time of 
960 s. The gray-scale shows the signal strength increasing from white to black (S/N range 2 ~ 20). 



2. Observations and data analysis 

During 1999, three observations of PSR J0030+0451 were 
performed with the 64m Parkes radiotelescope at a central 
frequency of 436 MHz. The total bandwidth was 32 MHz 
subdivided into 256 filter channels, each of width 125 kHz. 
The filterbanks are sampled continuously at a fixed rate of 
0.05 ms; the output from each channel is one-bit digitized 
and written to magnetic tape for off-line analysis. 

The first observation was performed on 1999 July 24 
and lasted 7000 s. The integrated pulsar signal was rela- 
tively strong, but unfortunately the observation time was 
short compared to the diffractive scintillation time and 
only part of the scintles are visible in the frequency-time 
dynamic spectrum. Since this MSP has a northern decli- 
nation, the longest observation time from Parkes is about 
5.5 hours, and that was the duration of two more obser- 
vations performed on October 15 and 16. In both cases 
the pulsar signal was weaker than in the first observation 
though we were able to detect a number of scintles in the 
dynamic spectra. 

During the analysis, pulse profiles are formed at the 
apparent pulsar period over short sections in time and for 
each of the 256 filter channels. Subintegration times were 
60 s for the first observation and 120 s for the other two. 
The signal-to-noise ratios in each subintegration for each 
frequency channel is then calculated. The first 16 channels 
were not considered in the analysis because of their mal- 
functioning; a few other channels that had bad response or 
in which strong interferences were present, were overwrit- 
ten by an average of the signal-to- noise ratios in adjacent 
channels. No smoothing operation was made on the dy- 
namic spectra but in all the observations a sum along the 
frequency and time axes was performed in order to im- 



prove the signal-to-noise ratio. Figure [l] shows a dynamic 
spectrum of PSR J0030+0451. 

To obtain the scintillation bandwidth and decorre- 
lation time-scale we performed a two-dimensional auto- 
correlation analysis on the dynamical spectra, following 
the method described by Cordes (1986). The scintillation 
bandwidth, A^, and the decorrelation time-scale, i<j, were 
obtained by fitting a Gaussian function to the zero lag in 
frequency and time respectively (details on the data anal- 
ysi s and descr iption of the software used can be found 
in Nigro 200C|) . Similar results were obtained with a 2- 
dimensional gaussian fit. 

3. Results and discussion 

Scintillati on velocities (Vj SS ) can be estimated using the 
equation (Gupta et al. 1994) 



(2) 



t; , or w 1A 4 V Av d DkpcX _j 

Viss = 3.85 x 10 — km s , 

while the line-of-sight-averaged electron density fluc- 
tuations can be estimated by flCordes et al. 1985 
|Cordes 1986; ) 



C% = 0.002 y^D-H / & A^ 6 m" 20 / 3 , 



(3) 



where Cjj is the amplitude of the (density fluctuations) 
power spectrum, Ai>d is the decorrelation bandwidth in 
MHz, td the characteristic time-scale in s, Dkpc the dis- 
tance in kpc, ^ghz the observation frequency in GHz and 
x = D /Dp is the scaling factor that accounts for an asym- 
metrically located scattering screen, whith D Q and D p the 
distances observer-screen and screen-pulsar, respectively. 
Apart (identifiable) peculiar cases, assuming x — 1 and a 
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Table 1. PSR J0030+0451 observation log and scintillation parameters. 





Date 


T 


V m 




tj 


V* 

* 1SS 








(hours) 


(km s _1 ) 


(MHz) 


00 


(km s _1 ) 


(m- 2 °/ 3 


24- 


-07-1999 


1.9 


11.7 


9.3 


10569 


12.5 ± 11.9 


-3.66 


15- 


-10-1999 


5.5 


29.2 


4.9 


10106 


9.4 ± 4.0 


-3.43 


16 


-10-1999 


5.5 


29.1 


2.0 


9518 


6.4 ± 1.6 


-3.10 



* Errors are 1<t level from equation 



thin screen model (Scheucr 1968) gives good agreements 
between V pm and Vi ss (see below). 

Table |l| lists the scintillation parameters of the three 
observations. Column 2 to 5 give observation time, trans- 
verse Earth velocity toward the pulsar, decorrelation 
bandwidth, characteristic time-scale, scintillation velocity 
and the measure of the turbulence along the line of sight. 
The statistical error on Vj ss is a combination of two quan- 
tities: the uncertainties on Avd and td from the gaussian 
fitting, and the empirical error of the scintillation mea- 
surements given by Bhat et al. (1998) as 



fd 



Ah> d t d 



-0.5 



(4) 



where a is the fractional error, B b s and t bs are the ob- 
servation bandwidth and time respectively, and fd is the 
filling fraction for number of scintles which here is as- 
sumed to be 0.5. This estimate uncertainty is far larger 
than that from the fit and that due to the dynamic spec- 
tra discretization, so these two latter are assumed to be 
negligible. 

The scintillation velocity is a combination of (pro- 
jected on the plane of the sky perpendicular to the di- 
rection of the pulsar) the pulsar velocity V pnl , ISM ir- 
regularities velocity ay ISM (assumed to be ~ 10 km s" 1 ) 
and Earth orbital velocity V® (varying from to ~ 30 



km 



In the case of a binary system, the pulsar's or- 



bital velocity should also be considered. For pulsars with 
high values of Vi ss , ay ISM and Vq give a negligible con- 
trib ution to the estim ate of the transverse velocity (see 
e.g. |Gupta et al. 1994|). For PSR J0030+0451 the values 



of Vi SS are similar to ov ISM and in two, out of three cases, 
it is Viss < V®. Therefore these velocities contribute sig- 
nificantly to the uncertainty. A possible way to consider 
these systematics is by 



ay, 



°Vis, 



Vis 



VA 



(5) 



However, to better estimate the velocity V± of the pulsar 
and the related uncertainty, a more reliable way would be 
performing a campaign of measurements preferably when 
the contribution of V® is minimum (~ 5 km s" 1 ). In ad- 
dition, having several measurements would allow to as- 
sume Vism = 0. We have planned to do so. Neverthless 
our measurements clearly indicate PSR J0030+0451 is a 
slowly moving object. 



3.1. LISM 

The average amplitude of the electron density power spec- 
trum value log = -3.33 we found for PSR J0030+0451 
is somewhat greater than —3.5, the average value for the 
whole pulsar sample with distances <; l kpc, i.e. out of the 



LISM (Cordes 1986). This estimate adds to that of other 
~ 50 pulsars with DM < 35 pc cm" 3 (20 by Bhat et al. 
(1998) and 31 by Johnston et al. (1998)). Most of these 
pulsars are localized within the LISM and have high 
values, confirming the existence of local enhanced scatter- 
ing. The overall picture shows clearly that the LISM is 
highly turbulent compared to the outer regions. Bhat et 
al. (1998) developed a model which they claim can fit the 
C^, observed in their pulsar sample. It is interesting to 
check their model against PSR J0030+0451 results since 
it is lo cated in the quadran t where no other MSPs are 
known ( Toscano et al. 1999a ). 

The expected decorrelation bandwidth value, given by 



Av d .j 



ll4 2 H>kpc /5 = 7-2 MHz, 



(6) 



is not too different from the mean measured value Avd = 
5.6 MHz. We computed the so called "anomaly param- 
eter" (defined a s A dm = (i%/^d 2 )obs/(*W^d 2 )cxp, see 
|Bhat et al. 1998|) for PSR J0030+0451 with respect to 
PSR J1133+16 and PSR J1929+10, from the Bhat et 
al. (1998) sample, and PSR J1744-1134 and J2124-3358 
(both MSPs) from the Johnston et al. (1998) sample, as 
they all have similar distances and, more importantly, 
DMs within ~ 1 pc cm -3 . The low values of A^m we obtain 
(1.4, 2.4, 1.60 and 1.34 respectively) do not confirm the 
decay trend shown in Figure 7 of Bhat et al. (1998). An 
asymmetry in the distribution of electrons in the LISM (as 
suggested by Lbscano et al. 1999a ) or a non-Kolmogorov 
density spectrum could explain the discrepancy. This is 
also shown in Fig. || where we plot yld m vs DM for 35 
PSRs (16 MSPs) with DM < 35 pc cm" 3 . This sample 
includes 32 objects from Johnston et al. (1998) and the 
3 MSPs: PSR J0030+0451, J1012+5307 and B1257+12. 
The plot was made by averaging the A dm of pulsars with 
DM differences less than 3 pc cm -3 and increasing the 
DM by that of the pulsars instead of a fixed step of 1 



pc cm 3 (as is done by Bhat et al. 1998). Error bars are 



inversely proportional to the number of pulsar products 
averaged. The three very high peaks at DMs 10, 21 and 
32 pc cm" 3 are due to the large discrepancies between 
the expected decorrelation bandwidths and the measured 
ones for PSRs J1730-2304, J1825-0935 and J2129-5718, 
respectively. The estimators Bhat et al. (1998) use to de- 
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fine the best model for the local scattering structure are 
in our case e^(DM) = 0.119 and e& = 0.275, confirming 
the inhomogeneity of the LISM. 

However the high variability of the observed decorre- 
lation bandwidth in PSR J0030+0451 by a factor - 5 is 
rather unusual. This is not caused by refractive scintilla- 
tion because the dynamic spectra show no frequency drift. 
Instead it likely originates from anisotropies in the density 
or the magnetic field structure of the LISM. Further and 
more sensitive observations are required to put any strong 
constraint on the proposed model for the LISM by Bhat 
et al. (1998). 




15 20 
DM (pc cnT 3 ) 

Fig. 2. "Anomaly factor" vs DM for 35 PSRs with 
DM < 35 pc cm~ 3 . The 3 isolated highest points are due 
to PSRs J1730-2304, J1825-0935 and J2129-5718. 



3.2. Transverse velocity 

The scintillation parameters allow us to derive a reliable 
estimate of the transverse velocity of PSR J0030+0451. 
We will show that the correlation between Vi ss and Vp m 
is very high and a relation Vi SS = V^, m can in general be 



assumed (see also pupta 1995[ ). As already mentioned, in 
spite of an average high flux density measured during the 
first observation, because of its short duration compared 
to td 1 the Vi SS estimate is affected by a relatively large 
error. Further, the data of the second and third observa- 
tion have lower statistical significance, so they can provide 
only upper limits (see also Tab. [l]). This suggests to give 
similar weights to the three observations in the calcula- 
tion of M ss . Hence we have determined a mean velocity 
of V± ~ 9 ± 6 km s _1 (la error). This number makes 
PSR J0030+0451 the pulsar with the smallest transverse 
velocity so far measured. Its position in the V pm vs. V^ ss 
scatter plot is quite isolated in the bottom left corner well 
in agreement with a correlation V pm = Vi SS (see Fig. ||). 

If we assume that the P = (1±0.2) x 10 -20 reported by 
Lommen et al. (2000) is all due to the Shklovskii effect, 



we obtain an upper limit to the transverse speed V± < 
65 km s . Now assuming Vj_ = 9 km s , the derived 
contribution to the spin-down rate is P, < 2 x 10~ 22 ss" 1 . 
This means that the Shklovskii effect introduces an error 
of only ~ 2% on P = P t + P s , lower than the ~ 20% so 
far obtained by the timing data fit. Hence we can assume 
the measured spin-down rate is the intrinsic one: P, = 
(1.0 ± 0.2) x 10" 20 s s" 1 ([Lommen et al. 20000 - However, 
using a 3er upper limit V± ~ 27 km s _i , we obtain a 
P s comparable to the P error estimate. We then conclude 
that Pi = (0.8-f-l) x 10~ 20 s s _1 , giving a surface magnetic 
field B = (2 -r- 2.2) x 10 8 G and a characteristic age r = 
(8 -T- 10) x 10 9 years. 

This pulsar was also detected in the soft X-ray band 



L 



( pecker et al. 2000 ). Its X-ray luminosity (0.1-2.4 keV) is 
(1 4- 2) x 10 30 (d/0.23 kpc) 2 erg s _1 . Using the in- 
ferred rotational energy loss rate E — (2.7-4-3.4) x 10 33 erg 
s _1 we obtain a relatively narrow X-ray efficiency range 
r) = L x /E ~ (0.3 4- 0.7) x 10- 3 (d/0.23 kpc) 2 . This is 
i n good agreement with t he proposed law L x /E ~ 10 -3 
(Becker & Triimper 1997) and points toward a magne- 
tospheric origin for the high energy radiation, in spite 
the inferred rotational energy loss rate and pulse period 
would put PSR J0030+0451 in the class of MSPs for 
which the X-ray emission is believed to be mostly ther- 
mal ( |Kawai fc Saito 1999| ). Becker et al. (2000), from the 
observed X-ray pulse profile (very similar to the radio one) 
and the large pulsed fraction (~ 69 ± 18%) also suggest a 
non-thermal origin of the X-ray emission. Thanks to the 
Chandra and XMM-Newton X-ray satellites, a good qual- 
ity X-ray spectrum should soon be able to give a definitive 
answer in this respect. 

PSR J0030+0451 is the 8 th isolated MSP detected for 
which a transverse velocity is available. Using for PSR 
J0030+0451 V pm = 32 km s" 1 (i.e. ~ 65/2), the VpZ for 7 
MSPs (we exclude PSR J1730-2304) is 58 ±6 km s" 1 (for 
5 MSPs Viss — 45 km s~ x ). This, even if at low statistical 
significance, is nearly a factor 2 lower than the 96 ± 9 km 
s" 1 foun d for 18 binary MSP s (for 12 MSPs Ws - 113 
km s" 1 ) fljohnston et al. 199$ [Toscano et al. 1999b|). 



4- V pm -Vi SS correlation 

Gupta (1995) made a comparison between proper motion 
and scintillation velocity for 59 pulsars using the Taylor & 
Cordes (1993) model to calculate their distances and a new 
formula to calculate the ISS velocities (Gupta et al. 1994). 
The comparison showed a good correlation between V pm 
and Vi SS . The correlation coefficient he found is r = 0.42. 
It becomes 0.47 if pulsars for which only upper limits were 
available are excluded from the sample. A least square fit 
gave a slope of 1.1 ± 0.13. 

We updated the Vp m -Vi ss diagram 

1. adding the proper motion and scintillation speeds of 



16 MSPs ( 


Toscano et al. 1999b 




Lange et al. 2000 


Johnston et al. 1998; 


Gothoskar & Gupta 2000 
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Fig. 3. Scatter diagram for proper motion and scintillation speeds for 77 pulsars. The grey points are the pulsars 
with updated values compared to Gupta (1995). The V pra = V[ ss correlation line is shown. 



Lorimer et al. 1995. this paper) and of the relativistic 

pulsar B 1534+1 2 flBtairs et al. 1998j ), 

using 14 revised proper motion estimates 



flFomalont et al. 1997 
using revised values 



Fomalont et al. 199S) 



of K s 



for 



(ghat et al. 1999j; |Lumiclla 2000 ) 



22 pulsars 
plus PSR 



B1259-63 flMcClure-Griffiths et al. 1998[ ). 



The V pm vs. Vi ss values for 77 pulsars are shown in the 
scatter plot of Fig. || and listed in Tab. ^, (MSPs are listed 
at the bottom of the Table). The uncertainties on Vi SS were 
computed according to Bhat el al. (1998) except for PSR 
B1257+12 and B1534+12 for which we gave a 20% er- 
ror. For the uncertainties on V pm we took into account 
the proper motion errors only. We know that uncertain- 
ties in the distance from the Taylor & Cordes model are 
on average ~ 30% but we do not include this source of er- 



ror. Only in a few cases do distances come from parallax 
measurements. 

For PSR J0030+0451 we adopted a (low statistical sig- 
nificance) proper motion value of V^m = 32 ± 32 km s _1 
(from the Shklovskii upper limit) and VJ SS = 9 ± 6 km 
s _1 (that is the average of the values of Tab. 0). The ob- 
ject at the top of the plot is PSR J2225+6535. Toscano 
et al. (1999a) gave for PSR J1730-2304 a lower limit of 
V pm iS 51 km s _1 and an upper limit of V pm Ss 107 km 
s _1 from the Shklovskii effect. We averaged these value 
and used the estimate of 79 ±28 km s _1 . For the relativis- 
tic pulsar PSR B1534+12 we used the inferred distance 
D = 1.1 kpc obtained by Stairs et al. (1998) assuming the 
validity of general relativity theory. 

A correlation analysis gives for the total sample a co- 
efficient r — 75%, much higher than the value obtained by 
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Table 2. Proper motion and scintillation velocities for 77 pulsars. 



PSR 




SVpm 


v iss 


5V iss 


PSR 


Vpm 


SVpm 


V iBS 


8V iBB 


(J2000) 


(km s" 1 ) 


(%) 


(km s" 1 ) 


(%) 


(J2000) 


(km s" 1 ) 


(%) 


(km s" 1 ) 


(%) 


0304+1932 


167 


21 


189 


15 


1907+4002 


126 


27 


1179 


11 


0323+3944 


237 


23 


163 


14 


1913-0440 


< 372 




380 


14 


0332+5434 


145 


7 


186* 


9 


1921+2153 


130* 


27 


173* 


9 


0358+5413 


135* 


26 


488 


18 


1932+1059 


71* 


8 


159* 


12 


0452-1759 


320* 


62 


96* 


20 


1935+1616 


480* 


21 


418 


13 


0454+5543 


202 


12 


404 


21 


1939+2134 


< 146 




205 


25 


0528+2200 


228 


99 


375 


22 


1946+1805 


37 


71 


61* 


12 


0543+2329 


376 


47 


574 


14 


1954+2923 


87 


45 


63 


27 


0614+2229 


< 304 




37 


20 


1955+5059 


495 


12 


1175 


12 


0629+2415 


< 538 




78 


48 


2018+2839 


< 31 




56* 


16 


0630-2834 


390* 


38 


170* 


8 


2022+2854 


97 


23 


250* 


7 


0659+1414 


251* 


19 


176 


28 


2022+5154 


91* 


2 


239 


20 


0814+7429 


75 


18 


89 


32 


2046-0421 


< 532 




328 


12 


0820-1350 


376 


13 


295 


16 


2046+1540 


159 


54 


202 


12 


0823+0159 


< 127 




53 


23 


2048-1616 


355* 


4 


501* 


5 


0826+2637 


196 


3.3 


279* 


22 


2113+2754 


381 


6 


1105 


24 


0835-4510 


141 


4 


136* 


20 


2116+1414 


< 618 




171 


26 


0837+0610 


174 


11 


214* 


16 


2157+4017 


485 


8 


1316 


8 


0908-1739 


142 


32 


394 


31 


2219+4754 


375 


31 


402 


22 


0922+0638 


919 


72 


642* 


4 


2225+6535 


1686 


2.3 


6772 


18 


0953+0755 


21 


28 


33* 


20 


2305+3100 


661 


28 


140 


25 


1115+5030 


142 


8 


235 


18 


0030+04511 


32 


100 


9 


66 


1136+1551 


475 


1.6 


396* 


40 


0437-4715 


121 


0.2 


170 


15 


1239+2453 


303 


4.4 


271* 


17 


0711-68301 


78 


2.5 


81 


20 


1456-6843 


89 


3.4 


41* 


8 


1012+5307 


62 1 


0.5 


50 


40 


1509+5531 


913 


5 


451* 


8 


B1257+12 


284 


1 


225 2 


20 


1543-0620 


< 170* 




80* 


8 


1455-3330 


100 


54 


81 


18 


1543+0929 


144 


40 


122 


13 


B1534+12 


132 


8 


161 


20 


1607-0032 


< 66 




60* 


8 


1603-7202 


27 


7 


81 


28 


1645-0317 


660 


42 


562 


16 


1713+0747 


35 1 


3.2 


29 


21 


1709-1640 


< 70* 




132 


19 


1730-23041 


79 


35 


62 


18 


1752-2806 


45* 


88 


55 


20 


1744-1134f 


36 1 


1.8 


31 


19 


1807-0847 


85* 


82 


95* 


11 


1911-1114 


183 


54 


180 


23 


1820-0427 


270 


13 


250 


10 


2051-0827 


14 


57 


100 


22 


1823+0550 


< 238 




62 


32 


2124-33581 


53 


3.8 


40 


23 


1825-0935 


75* 


60 


79* 


43 


2129-5718 


56 


50 


130 


26 


1840+5640 


286 


12 


529 


12 


2145-0750 


58 1 


11 


51 


17 


1844+1454 


481 


25 


424 


11 


2317+1439 


80 


16 


100 


20 



* Updated velocity estimates respect to Gupta (1995). f Isolated MSP. 
1 Lange et al. (2000). 2 Gothoskar & Gupta (2000). 



Gupta (1995). Further statistical analysis, using a rank 
Spearman test, gives a correlation coefficient r s — 79% 
at an extremely significance level of ~ 1 x 10~ 17 (the 
lower the value the higher the confidence of the result 
(Press et al. 1992)). A least squares fit with a linear law 
gives a slope of 0.81 ± 0.21. Changing the axes of the scat- 
ter plot (therefore taking into account the V pm error val- 
ues) we obtain a slope of 0.71 ± 0.20. Excluding the 14 
points with upper limits on V pm , we obtain r = 78% and 
r s = 80% at a significance level of ~ 1 x 10~ 15 and a 
slope of 0.88 ± 0.24. Exchanging the axes of the plot the 
least squares fit gives an angular coefficient of 0.70 ± 0.25. 
So we may estimate the slope of a straight line that fits 
the points taking into account the errors in both axis as 



m ~ 0.90 ±0.21 for the whole sample, and m ~ 0.82 + 0.25 
excluding the upper limits. 

4.1. MSP correlation 

PSRs J1603-7202 and J2051-0827 are the MSPs for 
which Vi ss differs significantly from V pm . To have Vp m = 
Viss the scaling factor should be x ~ 3.0 and x ~ 
7.1, respectively. These values suggest the presence of 
a scattering medium near the neutron star. This re- 
sult is in agreement with the theory that the binary 
system containing PSR J1603— 7202 evolved through 
a phase with critical unstable mass-transfer and most 
likely hosts a CO-WD companion rat her than a He- 
WD companion ( Tauris fc Savonije 1999 ). Binary systems 
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like this are believed to be surrounded by an enve- 
lope that slowly expands giving rise to the mass trans- 
fer (ran den Heuvcl 1994). This envelope is thought to 
be the main cause of scintillation. For the eclipsing bi- 
nary pulsar PSR J2051— 0827, the scintillation parameters 
vary significantly with orbital phase (Stappers et al. 1996; 
Johnston et al. 1998). It is likely that the wind blown from 



the companion star is contributing to the scintillation pa- 
rameters. In these cases, because the screen is so close to 
the pulsar, the measured scintillation velocity is effectively 
the velocity of the screen (see e.g. equation 3 in Britton et 
al. 1998). In both cases, the wind could easily have a ve- 
locity of ^100 km s _1 . One possible test of this idea would 
be to perform VLBI observations of these two pulsars to 
resolve (or not) the scattering disk; such an observation 
would give a direct measurement of the screen distance. 

We performed a correlation analysis between V pm and 
V^ ss for the sub-sample of 16 MSPs plus PSR B1534+12. 
Figure ^ shows the scatter plot. A simple correlation anal- 
ysis gives a coefficient r — 60%. A rank Spearman test 
gives a correlation coefficient r s — 67% at a significance 
level of ~ 3 x 10~ 3 . A least squares fit with a linear 
law gives a slope of ~ 0.57. Exchanging the axis of the 
scatter plot the least squares fit gives a slope of ~ 0.47. 
Excluding the two above mentioned MSPs from the analy- 
sis, the correlation coefficient becomes r = 87% and r s = 
91% at a much higher significance level of ~ 3 x 10~ 6 . 
A least squares fit with a linear law gives a slope of 
~ 1.1. Exchanging the axis of the plot the slope becomes 
~ 0.5. So we estimate the slope of the linear fit to be 
m ~ 1.1 for the whole sample, and m ~ 1.6 excluding 
PSR J1603-7202 and J2051-0827. 



100 




V pm (km/s) 



Fig. 4. Scatter diagram for proper motion and scintil- 
lation speeds for 16 millisecond pulsars. The V pm = Vi SS 
correlation line is shown. 



Since MSPs have low speeds, a rigorous analysis should 
take properly into consideration systematical uncertain- 
ties due to the Earth orbital velocity and the ISM irregu- 
lar motion. In spite of that the resulting Vp m -V; ss correla- 
tion is excellent. This demonstrates that averaging several 
observations spread over an Earth orbit is (in general) suf- 
ficient to wash out systematics. 



5. Conclusions 

We presented observational and derived parameters for 
the diffractive scintillation of the recently discovered PSR 
J0030+0451. Our results for this source show that: 

— its decorr elation bandwidth is highly variable; 

— the Kolmogorov coefficient logC^, = —3.33 along its 
line of sight is somewhat higher than the average for 
pulsars out of the LISM where it is log Cfj = —3.5. The 
enhanced turbulence and strong changes in v$ support 
the findings of Bhat et al. (1998) about the LISM, 
but seems to not be in agreement with their proposed 
model for the electron distribution; 

— it is unlikely its V± exceeds 30 km s -1 ; 

— as a consequence of this, the contribution of V± to the 
observed spin-down rate through the Shklovskii effect 
is negligible and suggests the derived parameters B, r 
and E are not affected by spurious effects. In partic- 



ular E = (2.7 -T- 3.4) x \Q iA erg s -1 and the observed 
X-ray luminosity are fully compatible with an X-ray ef- 
ficiency factor T) = 1 x 10 -3 . However the low observed 
L x suggests this pulsar is somehow in between the 
two proposed classes of X-ray emitting pulsars (ther- 
mal polar caps and magnetospheric emission). Only 
high quality X-ray spectral data will allow to discrim- 
inate between these two emission mechanisms. The X- 
ray new generation observatories XMM-Newton and 
Chandra can provide these data. 

To obtain accurate estimates of the scintillation param- 
eters, very long observation times and large band widths 
are needed. Because of its high sensitivity and site loca- 
tion in the North, the Indian Giant Metre Radio Telescope 
(GMRT) is the most appropriate instrument for PSR 
J0030+0451 scintillation studies. For minimizing system- 
atical errors, the observations must preferably be carried 
out at epochs when the projected Earth transverse orbital 
speed is low. Within a few years we should also be able 
to have an accurate timing measurement of V pm , although 
our results already show that PSR J0030+0451 is one of 
the slowest pulsars (with known V±). 

The pulsar sample for which proper motion and scin- 
tillation speeds are available was revised with new values 
for ordinary pulsars and 17 (recycled) millisecond pulsars. 
We performed several different correlation analysis on the 
resulting sample of 77 objects and found: 

— assuming the very simple model with the scattering 
screen placed midway between the Earth and the pul- 



sar, the overall V pm -Vi ss correlation is 



80%; this 
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value is much higher than the one reported by Gupta 
(1995); 

— MSPs show a good correlation r s ~ 65% that be- 
comes ~ 88% if two anomolous PSRs (J1603-7202 and 
J2051— 0827), for which local enhanced electron densi- 
ties were reported, are excluded. This confirms that 
interstellar scintillation parameters can be used as a 
probe of enhanced electron densities along the line of 
sight to PSRs; 

— we confirm that multi-epoch scintillation observations 
wash out systematics and give reliable estimates of the 
pulsar transverse velocity. 

Acknowledgements. F.N. wish to thank IRA-CNR for tech- 
nical support. The Australia Telescope is funded by the 
Commonwealth of Australia for operation as a National 
Facility managed by the CSIRO. 
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